Z2 topological insulator of ultra cold atoms in 
bichromatic optical lattices 
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We investigate the effect of a strong bichromatic deformation to the Z2 topological insulator in 
ultracold atomic system proposed by B. Beri and N. R. Cooper, Phys. Rev. Lett. 107, 145301 (2011). 
Large insulating gap of this system allows for examination of strong perturbations. We conclude 
that the Z2 topological character of the system is robust against a large global perturbation which 
breaks the inversion symmetry but preserves the time-reversal symmetry. 
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I. INTRODUCTION 

Topological insulators arc bulk insulators that have 
metallic states at their boundaries P, Robustness 
of these states against disorder and perturbations makes 
them promising for applications such as spintronics Q 
and topological quantum computation Topological 
invariants of the bulk material are essential for the ro- 
bust boundary modes. This urged consideration of topo- 
logical insulators on different lattice geometries [51-Q. 
Testing these systems against strong disorder and per- 
turbations is challenging both experimentally and theo- 
retically. Detailed numerical investigations reveal that 
topological protection is safe only for surface disorder 
smaller than the bulk band gap |10 |. 

Two and three dimensional topological insulators with 
band gaps in the order of the recoil energy have recently 
been proposed in ultracold atomic gases [ll| . This large 



gap systems are described in the nearly free electron 
limit, which is a different perspective to examine topolog- 
ical insulators compared to earlier studies [H, . The 
proposal utilizes interactions which preserves time rever- 
sal symmetry (TRS), analogous to synthesized spin-orbit 



coupling [12[, so that the insulators are classified by the 



so called I12 topological invariant 

Our aim in this Brief Report is to examine robustness 
of two dimensional Z2 topological insulators [ll| against 
large global perturbations that breaks inversion symme- 
try but preserves TRS. Motivated by their recent use in 
quantum simulations of relativistic field theories [isl [l6j 
and disorder induced localization [T7l - [l9j . we specifically 
consider here bichromatic optical lattice perturbations. 
Similar studies of Z2 topological insulators in bichro- 
matic potentials, but for a one-dimensional optical lattice 
in tight-binding limit, reveals that topological properties 
of the system can be probed by density measurements 

MM- 

This report is organized as follows. In Sec.|ll]we briefly 
review the proposal of Z2 topological insulator in ultra- 
cold atomic gases and introduce the bichromatic 
deformation of the lattice potential terms. In Sec. IIIII 
we describe the method of calculation of the I12 invari- 



ant in the case of absence of inversion symmetry [22l . |23[ 
and present the corresponding results. Wc conclude in 
SecHVl 



II. MODEL SYSTEM 

The Hamiltonian of an atom with N internal states 
with position r and momentum p can be written as. 



(1) 



here ^(r) is a position dependent potential which is a 
N X N matrix acting on the internal states of the atom. 
A Z2 topological insulator is invariant under the action 



of time-reversal operator Q 



icTyK, 



where ax y z are 



Pauli matrices acting on electronic spin and K is complex 
conjugation operator. This requires N to be even. The 
smallest potential matrix has = 4 and can be written 



^^^>-- \ Cl2 + ia-D {A-B)l2 



(2) 



Here A,B and C are real numbers, D is a 3-vcctor with 
real components and I2 is 2 x 2 identity matrix. 

This Hamiltonian can be realized us ing an atom with 
four internal states as proposed by Ytterbium 
^^^Yb, which has nuclear spin / = 2, is a good candidate 
for this purpose. It has a 2— fold degenerate ground state 
(^5*0 = g) and a long-lived 2— fold degenerate excited 
state ('^Po = e); furthermore there exists a state depen- 



dent scalar potential for a specific wavelength An 



in which the potential changes sign and becomes ±yam(r) 
for the ground and excited states. 

All the e — g transitions, shown in Fig. [Tal have the 
same frequency luq ~ {Eq — Eg)/h. We take the elec- 
tric field of the lasers interacting with the atom to be 
E = ee~'"* + e*e*'^*. The potential in the rotating wave 
approximation then can be expressed as [26j , 



V 



(|A + V;m)l2 -ia.tdr 
ia.edr -(fA + Km)ll2 



(3) 



where A = w — wq and dr is the reduced dipole moment 
of the atom . 

For the two-dimensional triangular lattice, following 
forms are assumed for the elements of the potential ma- 
trix, 

dre = [y(S,Fcos(r.ki),l/cos(r.k2)], (4) 



-A + Km(r) = ycos(r.(ki+k2)) 



(5) 



with ki = fc(l,0,0) and ks = fc(cos(6'), sin(6'), 0). Since 
w ~ wc have k ~ 2tt/ Xq with Xq — 578 nm the wave- 
length of the e — g transition. The spatial dependence of 
Kim is set by a standing wave at the antimagic wave- 
length Aam which creates a state-dependent potential 
with |ki-fk2| = 47r/Aam that leads to 6* = 2 arccos(27r/3). 
For simplicity, in all the following discussions we fix 
9 ~ 27r/3 and define a = 4:n/{^/3k). Therefore the op- 
tical potential has the symmetry of a triangular lattice 
with lattice vectors ai = (-^3/2, — l/2)a and a2 = (0, l)a. 

In order to examine the robustness of the Z2 topo- 
logical insulator to strong perturbations that breaks the 
inversion symmetry in the lattice, we introduce bichro- 
matic deformations. For simplicity, we do not change the 
angular orientation of the lattice so that the angular vari- 
able is not affected by the deformation. This can be 
accomplished by not deforming the element introduced in 
Eq. O In addition, we do not deform the first component 
in the element Eq. 31 The S dependent component plays 
the role of the spin orbit coupling to hybridize the lower 
bands of the two optical flux lattices. We choose not 
to perturb this essential process that mixes these bands 
with the non-trivial Chern numbers. 

The bichromatic deformation is introduced in the re- 
maining two terms such that 

drc = [VS, V cos(r.ki) -I- V' cos(2r.ki + 0), (6) 



Kcos(r.k2) + V" cos(2r.k2 -f (/)')]■ 

Here we examine the simple case where V = V" and 
= 0' = 0. This model is unitarily equivalent to the 
superposition of two bichromatic optical flux lattices that 
is coupled with a spin-orbit type interaction fllj . Such 
a configuration is illustrated in the Fig. Ilbl which is the 
cross section view of the potential along the ki direction. 
The small bumps at the potential minima is a global 
perturbation which is as strong as the insulating gap (see 
Fig. (2) • We explore its effect on the topological character 
of the insulator in the next section. 



III. RESULTS 

We determine the topological character of the insula- 
tor directly using the band structure to calculate the I12 
invariant. The Hamiltonian in Eq. ([T|) can be numeri- 
cally diagonalized using plane-wave method. In Fig. [5] 
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FIG. 1: (Color online) (a) Optical transitions between the 
ground (g) and excited ( e) energy levels. Two-fold nuclear 
spin degeneracy yields the manifold of states |gt)i \gi)j |et) 
and |e4,) used to define the potential matrix in Eq.[3](b) The 
schematic view of the cross section of the bichromatic optical 
lattice along the ki direction, for the case V' = 0.2V . 



the band structure is shown. In the case of inversion 
symmetric materials, finding the Z2 can be done using 
the Fu and Kane method [22| . Since the parity and time- 
reversal operators commute, we only need to check the 
parity eigenvalues in four time-reversal points in fc-space 
for all bands under the gap. The result for the system, 
without bichromatic deformation, turns out to be non- 
trivial so that it is either a Z2 insulator or a quantum 
spin Hall insulator, depending on the presence or absence 
of 6 coupling, respectively 

We can not use here the parity based calculation of Z2 
invariant given in Ref. 



ll| due to lack of inversion sym- 
metry in our case of bichromatic optical lattice. A gen- 
eral set of expressions of the Z2 invariants are developed 
in Ref. [22j • Equivalent expression based upon the B erry 
phase and the Berry curvature is constructed in Ref. [23 . 
This allows for a calculus recipe to determine the Z2 in- 
variant for a two-dimensional noncentrosymmetric mate- 
rial using the bulk band structure. We closely follow the 
notation and description given in Ref. [l^] to apply the 
method of Ref. [23| to our case. 

Time reversal constraint on the Bloch functions, as- 
sociated with the Kramer's degeneracy, allows for their 
determination using only the half of the Brillouin zone; 
as the other half is fixed by the constraint. In addition 
the constraint leads to a nonzero Z2 invariant as an ob- 
struction to smoothly defining the Bloch functions over 
the half Brillouin zone. Ref. [22| introduced the integral 
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(b) 

FIG. 2: (Color online) (a) Lowest energy bands for the lattice 
with V = 0.5Er are shown for the cases of (a) V' = and (b) 
V = 0.2V. The k-points are labeled as (Fp, = (pki + gk2)/2) 
and indicated in the figure by (p, q) with p,q (z {0, 1}. At the 
special points (1,0), (1,1) and (0,1) the 4— fold degeneracy 
is reduced to 2— fold by the presence of J = 0.25i5i{. All 
quantities plotted are dimensionless. 
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FIG. 3: (Color online) The integer n(kj) field for a 10 x 10 
grid representing the torus map of the Brillouin zone, where 
the shaded region is 23"*", is illustrated for (a) V' = 0.091^ and 
(b) V' = O.IV. Empty (white) circles, filled (black) circles 
and the blank rectangles correspond to n = — 1, +1 and n = 0, 
respectively. We denote ki = k2 — k. 



form of the Z2 invariant as, 

= dk-A-f (fkF{k) mod 2, (7) 

27r Jd<B+ J<3+ 

where 03"*" is half of the Brillouin zone, and A(k) = 
j ^^(u„(k)| VkU„(k)) is the Berry connection and F ~ 
Vk X is the Berry cm-vatm-c. The sum is over the oc- 
cupied bands which is n = 1, 2. ..4 for our case, where the 
lowest two bands are two-fold degenerate. The obstruc- 
tion of the Bloch wave function over the Brillouin zone 
gives the Z2 = 1 for topological insulator, and I12 = 
normal insulator. 

In order to implement this formula numerically, a link 
variable Ujkj) = det(u„(kj)|-u,„(kj + fj.)) /J\f(kj) is in- 
troduced [23], where A/'(kj) = |det(-u„(kj)|w„i(kj + fJ-))\ 
here n is the unit vector on mesh. The discrete version of 
Berry connection and corresponding Berry curvature are 
given as A =Im lnC/^i(kj) and F =Im \nU^(kj)U,^{k.j + 
li)U~^{'kj + i')U~^(kj) respectively, where the com- 
plex logarithm is defined within its principal branch 
(— 7r,7r]. Since we have time-reversal symmetry in our 
system, Kramer's pair at time-reversal invariant points 
0, ki/2, k2/2 and(ki +k2)/2 for £„ = £„_i need to obey, 

|ii„(k,)) = e|7i„_i(kj)). (8) 



After numerically evaluating the Bloch functions at each 
grid point on the Brillouin zone, we use this extra gauge 
condition for the time-reversal constraint to determine 
the Bloch functions at the critical points. 

The discretized version of Eq. ^ can be written as the 
summation over the half Brillouin zone D ~ ^j,. n(fcj) 
where the integer n(kj) field is defined for each plaquette 
of ^+ by 

n(k, ) - ^ [( A, - A. A^) - F(k, )] . (9) 

Here, A^j represents the forward differencing operation. 

To calculate the Bloch wave functions we need to di- 
agonalize the Hamiltonian Eq. ([T]) with the potential el- 
ements of Eq. ini The band structure is shown in Fig. [21 
Four-fold degeneracy of the lowest band is removed by 
the 6 dependent spin-orbit like interaction and reduced to 
two-fold degeneracy except at a few points. Upper panel 
shows the case when there is no bichromatic deformation. 
In this case there are three points, (1, 1), (0, 1) and (1, 0), 
which are four-fold degenerate. Under the bichromatic 
perturbation only the point (1,1) can retain its four- fold 
degeneracy as shown in the lower panel. By changing the 
values of V' from zero to values as large as the insulating 
gap, we see that in special points in fc— space the degcn- 
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eracy is lifted eventually. The gap becomes as large as 
for V > 0.2Eji and remains preserved despite the 
strong perturbations V as large as the gap itself. The 
number of occupied bands that are relevant for the Z2 
invariant calculation is still four when the bichromatic 
lattice is introduced. 

The points T„in shown in the band diagram corre- 
sponds to the upper right quadrant of the Brillouin zone 
depicted in Fig. |3l We determine the Bloch functions for 
the lowest four bands at every mesh point in the shaded 
region, *B+, of the Brillouin zone. The Bloch functions 
for the critical points are updated according to the time- 
reversal gauge constraint in Eq. |8l This final set of Bloch 
functions are then used to determine the Berry curva- 
ture and the Berry connection to calculate the integer 
field for the each plaqucttc of *8+. The sum over the 
half Brillouin zone *B+ gives the I12 invariant D. It can 
be checked in Fig. [3] that the sum over the entire zone 
gives the first Chern number which is zero. We repeat 
the calculation systematically for different strengths of 
the perturbation potential. The results arc summarized 
in Tabic n 

TABLE I: The D parameter for different values of strength of 
the bichromatic optical lattice's depth relative to the primary 
lattice. The and 1 values of the D parameter show the 
triviality and non-triviality of the system respectively. 



Vi/V 


D 
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0.01 





0.02 


1 


0.03 


1 


0.04 


1 


0.05 
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0.06 


1 


0.07 


1 


0.08 





0.09 





0.1 


1 


0.2 


1 


0.3 


1 



Table U indicates that the system holds its nontriv- 
ial character for majority of different V' cases. The oc- 
currence of few zeros for the topological parameter D is 
due to the numerical difficulties addressed in other places 
as well [2^. The difficulties in implementation of the 
present methods, especially for first-principles electronic 
structure calculations in two and three dimensions, moti- 
vate the search for gauge independent and easier methods 
to calculate Z2 invariants [231 ■ ^'^^ our case the method 
is relatively simple to use, and hence we did not use such 
alternatives. 

IV. CONCLUSION 



We consider the Z2 topological insulator of ultracold 
atomic g^tem on a triangular optical lattice, introduced 
in Ref. deformed globally by bichromatic potentials. 
We examine the effect of this deformation which can be 
as large as the insulating gap, on the topological char- 
acter of the system. We find that the system retains its 
topological character robustly against this kind of per- 
turbation. The large insulating gap of such topological 
bichromatic insulators with deep bichromatic deforma- 
tion can be used for disorder and relativistic dynamics 
studies in the nearly free electron regime. 
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